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Abstract 

Highly charged relativistic heavy ions have high cross-sections for two-photon interactions. The 
photon flux is high enough that two-photon interactions may be accompanied by additional pho- 
tonuclear interactions. Except for the shared impact parameter, these interactions are independent. 
Additional interactions like mutual Coulomb excitation are of experimental interest, since the neu- 
trons from the nuclear dissociation provide a simple, relatively unbiased trigger. 

We calculate the cross sections, rapidity, mass and transverse momentum (pr) distributions for 
exclusive 77 production of mesons and lepton pairs, and for 77 reactions accompanied by mutual 
Coulomb dissociation. The cross-sections for 77 interactions accompanied by multiple neutron 
emission (XnXn) and single neutron emission (lralra) are about 1/10 and 1/100 of that for the 
unaccompanied 77 interactions. We discuss the accuracy with which these cross-sections may 
be calculated. The typical pt of 77 final states is several times smaller than for comparable 
coherent photonuclear interactions, so pt may be an effective tool for separating the two classes of 
interactions. 

PACS numbers: PACS Numbers: 25.20.-x, 12.40.Vv, 13.60.Le 
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I. INTRODUCTION 



With their large charges Z, relativistic heavy ions carry strong electromagnetic fields 
which act as strong sources of nearly-real photons. These photons can induce a wide variety 
of photonuclear and two-photon physics [1]. The photon flux scales as Z 2 , so the two-photon 
cross-section scales as Z 4 . The particle production rates are competitive with those obtained 
at e + e~ colliders. At the LHC, there is strong interest in using two-photon physics to search 
for signs of new physics [2J. 

Because the coupling constant Za w 0.6 (a ~ 1/137 is the electromagnetic coupling 
constant) is large, two-photon interactions may be accompanied by additional electromag- 
netic interactions, such as additional two-photon reactions, or photonuclear interactions. 
Photonuclear reactions can produce collective nuclear excitations, e.g. a Giant Dipole Res- 
onance (GDR), or can involve more energetic processes, such as pion production. Multiple 
additional reactions also occur. For example, one photon from each nucleus can dissociate 

in 

the other, in mutual Coulomb dissociation (MCD) [3|, |4{]. Experimentally, neutrons from the 
nuclear dissociation are detected in zero degree calorimeters (ZDCs) downstream from the 
collision points, in both directions. These neutrons make a convenient trigger, separating 
two-photon and other ultra-peripheral collisions (UPCs) from backgrounds such as cosmic 
rays and beam-gas interactions. The neutrons also 'tag' events with a smaller mean impact 
parameter than exclusive two-photon events. Neutron tagging has been used to study e + e~ 
5, 6] and vector meson photoproduction ls[ ^]. 

Tagging is of special interest for e + e~ because the stronger fields at small impact pa- 
rameters, where non-perturbative or strong-field electrodynamic effects will be strongest. 
Calculations of e + e~ pair production accompanied by MCD are challenging, because the 
ypical impact parameters are smaller than the electron Compton wavelength, Ac* = 386 fm 



10j . The cross-section measured by STAR is larger than lowest-order QED calculation; the 



difference can be explained by higher order corrections [11] . 

Here, we calculate the cross sections, rapidity, invariant mass, and transverse momentum 
distributions for two-photon production of a selection of scalar and tensor mesons and 
and t + t~ pairs. We calculate the cross-sections and distributions with and without nuclear 
breakup, and explore how nuclear breakup affects the reactions. 

Because the cross-section for n-photon interactions scales as Z 2n we focus on heavy nuclei: 
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FIG. 1: The dominant Feynman diagrams for two photon reactions (a) without and (b) with nuclear 
excitation. Cross diagrams and, for (b) time reversed diagrams are omitted; due to factorization, 
they simply add. 

gold-gold collisions at a center of mass energy of sjs nn = 200 GeV per nucleon (beam Lorentz 
boost 7 = 108), as seen at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory, and lead-lead collisions at an energy of 5.5 TeV per nucleon, as will be 
produced at the large hadron collider (LHC) at CERN, where 7 = 2940. We assume that 



121. The 



RHIC runs for 10 s/year, at an average gold-gold luminosity of 2 x 10 /cm /s 
LHC is expected to devote ~ 1 month, or 10 6 s/year to lead beams, at an average luminosity 
of 10 27 /cm 2 /s Q. 

Section II will present a calculation of the two-photon luminosity. Section III discusses 
the effects of multiple interactions, while Section IV finds the two-photon luminosity under 
various different nuclear breakup conditions and discuss the uncertainties in the calculation. 
The final state pt spectra are discussed in Section V, and Section VI gives some conclusions. 



II. TWO-PHOTON LUMINOSITY 



According to the method of equivalent photons, the cross section for a two-photon reac- 
tion, Fi g, l a, factorizes into an elementary cross section for 77 — > X and a 77 luminosity 



C 



77 



13, 14j. The cross section to produce a final state with mass W is 

n(ki) n(k 2 ) 



a(A+ A^A+A+X) 



dk\dk<i 



-(7(77 -> X(W)) 



(1) 



h k 2 

where k± and k 2 are the two photon energies, and n(k) is the photon flux at energy k. 
The 77 luminosity is given by convolution of the equivalent photon spectra from the 
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two nuclei. In impact parameter space, the total number of photons from one nucleus is 
obtained by integrating over all impact parameters larger than some minimum, typically 
given by the nuclear radius. This is similar to integrating over all possible momentum 
transfers, Q, from some minimum determined by the kinematics up to a maximum given 
by the nuclear form factor. In hadronic collisions, the impact parameter representation 
provides the best way to incorporate effects of strong absorption. Hadronic interactions will 
dominate in collisions where both hadronic and electromagnetic interactions are possible. 
The hadronic interaction must therefore be excluded to obtain the effective or usable cross 
section for two-photon interactions. In impact parameter space, this can be accomplished 
by requiring that the impact parameter be larger than the sum of the nuclear radii. The 



equivalent two-photon luminosity is thus 



m 



dL - -£aaZ[ Al I d 2 b 2 n(k 1 ,b 1 )n(k 2 ,b 2 )Q(\b 1 -b 2 \-2R A ) (2) 

Z J b\ >Ra Jb 2 >R A 

where Caa is the ion-ion luminosity, n(k, b) is the flux of photons with energy k at impact 
parameter b, and Ra is the Woods-Saxon nuclear radius. 

The requirements b\ > Ra and b 2 > Ra ensure that the final state is produced outside 
the nuclei. Otherwise, the final state will usually interact with the nucleus, destroying itself 
and breaking up the nucleus. This requirement may not be strictly necessary for some final 
states, such as lepton pairs. Alternately, a smaller radius might be appropriate. However, 
since the electric fields drop rapidly for b < Ra, this is a relatively small correction. The 
function imposes a requirement that the nuclei not physically collide; we discuss more 
detailed hadronic interaction models below. 

The photon flux may be modelled using the Weizsacker- Williams method. For 7 ^> 1 

d 3 N Z 2 a . . . 

where x = bk/j. Here, K\ is a modified Bessel function. The two photon energies k\ and 
k 2 determine the center of mass energy W and rapidity y: 

W 

h,2 = ~e ±y (4) 

and 

y = l/2Hh/k 2 ). (5) 



dWdy 
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The maximum effective two-photon energy, W max occurs at y = 0, when k\ = k 2 = j/Ra- 
W m ax is about 6 GeV for gold at RHIC, and 150 GeV for lead at the LHC. W max is higher 
for lighter nuclei and protons, but the 77 luminosity per collision is lower, and multi-photon 
interactions are very rare. 

Equation (j2J) treats the nuclei as hard spheres with radius R A - Since there is a finite 
probability to have hadronic interactions at impact parameters b > 2Ra, more accurate 
hadronic interaction probabilities can be included by modifying Eq. (j2J), to 

£ 77 „ W 



dWdy 



Caa^t I d\ [ dXnihMMk^b^l-PHdh-hl)} (6) 

Z Jb,>R A Jbn>R A 



with the hadronic interaction probability 



P H (b) = l-exp^-a nn J d 2 f T A (f) T A {r- S) J . (7) 

a„„ is the total hadronic interaction cross section, 52 mb at RHIC and 88 mb at the LHC 
12|. We use the total cross sections, since even an elastic nucleon-nucleon interaction will 
break up the nucleus. The nuclear thickness function is the integral of the nuclear density, 
p(r) 

T A (b) = / dzp(b,z)dz (8) 



where b is the impact parameter from the center of the nucleus. The nuclear density profile 



p(r = y\b\ 2 + z 2 ) of heavy nuclei is well described with a Woods-Saxon distribution. We 
use parameters determined from electron scattering data (R=6.38 fm for Au and R=6.62 
fmforPb) 3- 

As Fig. [2] shows, Eq. ([6]) gives 77 luminosities about 5% lower than the hard sphere 
model for W = 0.1W max , falling to 15% lower for W = W max [n] ]. These differences are 



comparable to those found elsewhere [16 ]. 



Baur and Ferrara-Filho derived Eq. (j2J), and then used a change of variables and the 
function to reduce the dimensions of the integral [l?j]. Although this approach speeds the 
calculation, it works poorly for the realistic models of nuclear density or when additional 
photon exchange is included. Cahn and Jackson used a related approach, calculating the 
luminosity analytically without the requirement b > 2Ra, and then numerically calculating 
a correction for the overlap [3]. This approach also cannot accomodate nuclear breakup. 

In principle, the meson or lepton pair production cross section depends on the angle 
between b\ and b 2 , with different cross sections for parallel and perpendicular photon po- 
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FIG. 2: Reduction in two-photon luminosity when the nucleon density is approximated with a 
Woods-Saxon distribution (dotted line) as compared with a flat distribution for hard sphere nuclei, 
requiring \b\ — &2I > (solid line). The calculation is for gold-gold interactions at RHIC. 



larizations [19(. However, after integration over b\ and 62, the differences are small, and we 
neglect them. 

The cross section to produce a narrow resonance with spin J and mass m is 
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87r\2J + l)^- 2 5(W-m 



(9) 



where we neglect the width of the hadronic resonances. The magnitude of the error due 
to this narrow-resonance approximation scales linearly with the width. For coherent pho- 
tonuclear p° production, including the width reduces the cross section by about 5% [l2j. A 
similar correction is expected for two-photon production; these adjustments should should 
scale linearly with the relative width (T/M) of the resonance. 

We also consider the production of continuum lepton pairs. For leptons with mass M, 
the cross section is given by the Breit- Wheeler formula 

W + VW 2 - AM 2 
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\M 1 f l 1M 2V1 



(10) 



W 2 W A J 2M V W 2 V W 2 

This equivalent-photon approach is simpler than a full QED calculation in that it neglects 
the virtuality of the intermediate photon lines in Fig. 1. For this reason, we do not consider 
e + e~ pairs, where this intermediate state affects the pair p T distribution jol. flo|. 



III. TWO-PHOTON REACTIONS ACCOMPANIED BY NUCLEAR BREAKUP 



Because Za is large, two-photon reactions may be accompanied by additional photonu- 
clear reactions, as in Fig. (lb). As long as the two meson-producing photons do not excite 
the nucleus that emitted them, this is the dominant diagram for producing a meson while 
exciting both nuclei. As with photoproduction of vector mesons accompanied by MCD 
the kinematics of photon emission does not favor nuclear excitation. Then, the individual 
photon reactions are independent 20j and the processes 

factorize fl. Photonuclear breakup 
can be incorporated in Eq. ([6]) by including the photonuclear excitation probability, P(b): 

£77 _ r W f a2 



dWdy 



Caa^t ! d\ I rf 2 6 2 n(fci,&i)n(A; 2 ,6 2 )P(6)[l-P^(6)], (11) 

Z J b\ >Ra Jb 2 >R A 



where b = \bi — 62 1- Here, P{b) is the total breakup condition of interest, whether it is 
one nucleus or two. The two-photon luminosity accompanied by single nuclear excitation 
was calculated in this approach in Ref. 22|. We separately consider two cases. The first, 
labelled Xn covers all nuclear excitation, including high-energy excitations which may emit 
pions in addition to dissociating the nucleus. The second, single neutron emission, In, is a 
subset of Xn, in which a single neutron is emitted. 

The lowest order probability for Coulomb breakup of a specific nucleus is 

p *^=L dk ^k a ^ k) - (i2 » 

where k is the photon energy and d 3 n^/d 2 bdk the photon density from Eq. (3). ct 7 a^a*(&), 
;he excitation cross section, is determined by data collected over a wide range of energies 



23| ; E min is the minimum energy for this excitation, while E max = ^Tic/b is the maximum 



rows that this is the 



24]. 



photon energy for which there is significant flux. The superscript 1 s 
lowest order probability. A similar calculation was performed in Ref. 

More precisely, Px n (b) is the mean number of excitations; the probability of having 
exactly N excitations follows a Poisson distribution, so the probability for at least one 
Coulomb excitation is 

P Xn (b) = l-ex P (-P 1 Xn (b)). (13) 

Single neutron emission usually occurs when the nucleus is excited into a giant dipole 
resonance (GDR), where the protons and neutrons oscillate against each other 25|. Photons 
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with energies from 8-24 MeV can excite a GDR resonance in heavy nuclei. To observe the 
In breakup, the GDR excitation must not be accompanied by any higher-energy reactions: 



P ln (b) = P} n (b)e W (-P 1 Xn (b)). 



(14) 



The lowest order probability Pl n (b) is determined as in Eq. (|12p . except that the photon 
energy integration is truncated at the maximum GDR energy, 24 MeV for both nuclei. Data 
on photoemission of single neutrons is used, avoiding uncertainties in the GDR branching 
ratios [26] . 



In MCD, the two nuclear breakups occur independently [4j, ]2l|, |23j , so the probability for 
MCD is 

PxnXn(b) = (Pxn(b)) 2 (15) 

and 

iW&) = (Pm(b)) 2 . (16) 

This independence has been verified by comparing the cross-sections for single nuclear 
breakup, both In and Xn, and MCD |27j]. It also appears to hold for p° photoproduction 
with accompanying MCD 



IV. CROSS SECTIONS, RAPIDITY AND MASS DISTRIBUTIONS 



Table I shows the two-photon widths and cross sections for the different tags at RHIC and 
the LHC, for several meson resonances and for and r + r~, from Eq. (TTTT) . For a 'stan- 

dard' RHIC year, the 550 fib cross-section for the / 2 (1270) leads to 1,100,000 events/year. 
At the LHC, a 1/zb cross section corresponds to 1,000 events/year. 
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Untagged meson production rates have been calculated by many authors 
Unfortunately, some authors used different ion species and Lorentz boosts. Several earlier 
papers considered uranium beams at RHIC with 7 = 100. These papers found rf production 
cross-sections of 2.9 mb 



35( to 3.6 mb [30] . For lead at the LHC, with 7 = 4000 (rather 
than the planned 7 = 2940 which is used here), they find n' cross-sections of 30 and 32 mb 



respectively, compared to the 22 mb found here. Roldao and Natale [3l| consider uranium 
beams at RHIC, and lead beams with 7 = 3360 at the LHC, using the Cahn- Jackson 
approach [181 ]; our RHIC results are about 40% lower than theirs; the LHC results are quite 
comparable, except that our r\ c cross-section is 6 times larger than theirs. It is difficult 



8 



TABLE I: Spin/parity, two-photon widths T 77 |28f] , and cross sections for two-photon production 
of various final states with and without nuclear breakup. 



Meson 




1 77 


overall 


XnXn 


lnln 






(keV) 


a [nib] 


a [jib] 


a [fib] 


Gold beams at RHIC 


n 


o-+ 


0.510±0.026 


1.05 


51.7 


4.9 


rf 


o-+ 


4.30±0.15 


0.72 


49.7 


4.6 


/ 2 (1270) 


2++ 


2.60±0.24 


0.55 


44.9 


4.1 


$(1525) 


2++ 


0.081±0.009 


0.0069 


0.63 


0.057 


T]e 


o-+ 


7.2±0.7±2.0 


0.0029 


0.40 


0.034 








142 


5,628 


537 


T + T - 






0.00079 


0.14 


0.011 


Lead beams at LHC 


n 


o-+ 


0.510±0.026 


18.8 


337 


21.6 


i 


o-+ 


4.30±0.15 


21.9 


469 


30.3 


/ 2 (1270) 


2++ 


2.60±0.24 


23.4 


562 


35.7 


$(1525) 


2++ 


0.081±0.009 


0.38 


9.7 


0.62 


Tie 


o-+ 


7.2±0.7±2.0 


0.57 


19.3 


1.2 








2,017 


33,084 


2,128 


T + T - 






0.55 


22.8 


1.4 



to understand the large difference for the r\ c . The agreement between these calculations is 
generally good, but detailed comparisons are not possible. 

Newer calculations have used the same species and boosts. Our cross-sections are to 

n 

20% higher than those in Ref. 36]; this work seems to use 7 = 100 for RHIC. The exception 
is r] c production at RHIC, where the current cross-section is 6 times higher. The difference 
may lie in how the papers handle the non-collision requirement; Bertulani and Navarra apply 
a form factor to the nuclei instead of applying an explicit requirement that b > 2R; this can 
lead to substantially different results near the kinematic limit. 

n 

These rates are also very similar to the r\ c results reported in Ref. [29]. Our rates are 
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10-25% higher than those in a very similar calculation [151 ]; the difference stems from slightly 
different Ra (Ref. [15] used an A-dependent parameterization of the nuclear radius, rather 
than electron scattering data [12|) and older values for the two-photon widths, T 77 . 

The RHIC fi + fi~ and r + r~ rates are about 20% and 33% respectively lower than in Ref. 



29( | . The t pairs are close to threshold, so are very sensitive to the detailed of the nuclear 



surface calculation. For the /i pairs, the difference may be because Ref. [29( integrates 
over rapidity rather than the impact parameter for the second photons; the effect on the 
nuclear and final-state overlap requirements are not transparent. The RHIC and LHC fi + fi~ 
rates are also about 37% and 22% respectively lower than approximate Born calculations 



08 and 3000). One 
, and found results 



of Hencken, Kuraev and Serbo [37j] at roughly the same energies (7 of 
of us has also performed Born calculations by numerical integration |3£ 
consistent with Hencken, Kuraev and Serbo when adjusted for the slightly lower energies of 
the numerical calculations (7 of 100 and 2760). These Born calculations have no cutoff in 
integration over h\ and 6 2 ; point charges are assumed for the ions and a correction is made 
by applying a form factor for the virtual photons emanating from each ion. 

Notable differences between the various calculations are in the treatment of the possibility 
of accompanying hadronic interactions and in how the virtual photons are cut off, either by 
integration over b\ and 62 or by a form factor. Here, the ultimate accuracy depends on a 
knowledge of the nuclear matter distribution; although there are accurate measurements of 
the proton distributions in heavy nuclei. Data on neutron radii is sparse, but it appears 
that neutrons are less confined than protons in heavy nuclei. For lead the neutron radius 
has been estimated to be 0.17 fm larger 32[. Figure [3] shows how £ 77 varies with W for 
a +/-10% change in nuclear radius, for gold beams at RHIC. A ± 10% change in radius 
leads to a ± 10-30% change in luminosity. Uncertainty in the nuclear radius leads to an 
irreducible systematic uncertainty in 77 luminosity. 

A similar limitation holds for 77 interactions with protons, where the choice of proton 
form factor or radius can substantially affect the photon flux 39f| . Efforts to use lepton 



pair production to measure the luminosity at the LHC 40j have to take into consideration 



these radius and form factor uncertainties as well as the possibility of higher order Coulomb 
corrections arising from the large Z heavy ions [2J. 

The XnXn and lnln cross-sections in Table I are new. They are about 1/10 and 1/100 
of the untagged cross sections, respectively, at RHIC; at the LHC the reduction is a bit 
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FIG. 3: Two-photon luminosity £ 77 as a function of W for two different nuclear radii (+/-10% 
change) for gold beams at RHIC. 

larger. These are small fractions, but XnXn reactions may still be of value because of the 
triggering advantages that they provide. 

The cross sections and corresponding production rates for single mesons are high. For 
example, one can expect to produce ~ 10 6 and ~ 23 • 10 6 /2(1270) mesons in a standard 
RHIC and LHC year, respectively. If one requires mutual Coulomb breakup (XnXn) the 
rates are reduced to ~ 90, 000 and ~ 560, 000. The n c might be hard to detect at RHIC, 
but the rates should be sufficient at the LHC, where around 600, 000 n c mesons should be 
produced in one year. Of these, ~ 20, 000 should remain if one requires XnXn breakup. 
Although these rates are high, it should be noted that the cross section for coherent vector 
meson production with similar event topology are about two order of magnitudes larger for 
mesons of similar mass 0, [ijj]. A way to separate the two reactions may be through the 
different pr-spectra, as will be discussed in the next section. 

Figure H] shows the rapidity distribution da/dy for /2(1270) and r] c production at RHIC 
and the LHC. Spectra for XnXn and lnln breakup are shown, along with the untagged 
da/dy. The distributions with breakup are slightly narrower than those without. This is 
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FIG. 4: The rapidity spectrum da/dy for /2(1270) with (a) gold beams at RHIC and (b) lead 
beams at the LHC, and for the rj c with (c) gold at RHIC and (d) lead at the LHC. Three curves 
are shown in each panel: the total two-photon cross section (solid curve), and with XnXn (dashed 
curve) and lnln (dotted curve) excitation. 

because MCD selects events with smaller impact parameters. This increases the fraction of 
events with hard photons, eliminating interactions at large impact parameter, where a low 
energy photon far from it's source nucleus interacts with a higher energy photon near it's 
source; the large energy disparity corresponds to a large rapidity. 

Figure [5] compares the fi + fi~ and t + t~ pair mass distributions at RHIC and the LHC. 
Subject to the aforementioned caveats about beam species and energies, the untagged spectra 
are similar to those found by other authors Q, 0, 41]. The solid line shows the untagged 



spectra, while the dashed and dotted lines are for XnXn and lnln, respectively. The mass 
spectra with breakup are harder than those without. 
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FIG. 5: The invariant mass spectrum da/dM for n + n~ pairs with (a) gold beams at RHIC and 
(b) lead beams at the LHC, and for the r + r _ with (c) gold at RHIC and (d) lead at the LHC. 
Three curves are shown in each panel: the total two-photon cross section (solid curve), and with 
XnXn (dashed curve) and lnln (dotted curve) excitation. 

V. TRANSVERSE MOMENTUM SPECTRA 

Ultra-peripheral interactions (UPCs) are fully coherent with both nuclei, and so have a 
very small final state pt', this characteristic is important for separating UPCs from back- 
ground events. Here we consider a more complicated problem, that of separating fully re- 
constructed 77 interactions from coherent photonuclear interactions. The final state meson 
Pt is the vector sum of the two photon perpendicular momenta, k±. We assume that, inte- 
grated over all of transverse space this angular distribution is flat, and add the two-photon 
k± in quadrature. 
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At a distance b from the center of the emitting nucleus the photon flux is [1 



d 3 n{b, k) Z 2 a 



(17) 



dkd?b TT 2 k 

where J\ is a Bessel function. 

For heavy nuclei, the nuclear charge form factor F can be analytically modelled by the 
convolution of a hard sphere with a Yukawa potential of range a = 0.7 fm 12] : 



4irp 



F (q) = -7T [sin(qR ws ) - qR ws cos(qR ws )] (— — — ) 
Aq d L 1 + a z q z 



+ a 2 q 2 ' 



(18) 



Here, p° is the density and Rws is the Woods-Saxon radius of the nucleus. This gives 
an excellent approximation to the Woods-Saxon distribution. The radii of the heavy ions 
are well measured however recently published results indicate that the neutron and proton 
distributions differ in the nuclei which may lead to a larger radii of the matter distribution 
and essentially limit the accuracy of the photon flux determination 32]. 

The k± spectra of the virtual photon fields are complicated because k± and transverse 
position are conjugate variables, and cannot both be defined at the same time. This compli- 
cates any determination of the k± spectrum with constraints on transverse position, such as 
those imposed by photoexcitation reactions. We avoid this problem by selecting the photon 
energies, and then determining the transverse momenta solely by using the energies, using 
Eq. (jTTI) . The photon k± spectrum for fixed k is given by [331 ]. 

dN _ 2Z 2 aF 2 (k\ + k 2 h 2 )kl 

dk^ ~ n[kl + k 2 / 1 2 } 2 ' ^ ' 

This is the same approach used to calculate px spectra for vector meson photoproduction 



331 ] . Hencken and collaborators used Eq. (|T9|) . with a point-particle form factor and a cutoff 



Pt < 1/Ra (221] • Other groups (including some of the present authorship) used this equation 
with a Gaussian form factor 34 ] . 

The photon k± distribution is asymmetric, with a maximum at k± = \f?>k/^ and large 
tail at high kj_. The distribution has two different momentum scales, %/Ra and k/j. For 
h/a > kj_ ^> Ti/Ra and k± ^> k/j, the tail falls as l/kj_, regulated by the form factor. 

The final state pt is the vector sum of the two photon kj_. We assume that the angle be- 
tween the two photons is completely random, neglecting the possibly different cross sections 
for photons with parallel and perpendicular spin. 
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FIG. 6: The transverse momentum spectrum dN/dpx at mid-rapidity, y = 0, for two-photon 
production of states with masses of 0.2 (solid), 0.5 (dashed), 1.0 (dotted) and 3.0 (dash-dotted) 



GeV in Au + Au interactions at RHIC. The curves are normalized to have equal areas 



42] 



Figure [6] shows the pt spectra at mid-rapidity, y = 0, for two-photon production of states 
with masses of 0.2, 0.5, 1.0 and 3.0 GeV in AuAu collisions at RHIC. The spectra peak 
around y/l.BM/'y. The average px is around 75% larger, because of the long high-p^ tail. 

In photonuclear interactions, most of the pt comes from the nucleus; the typical pr is is a 
few Ti/Ra, or 50-100 MeV/c jj, 9]. This is several times higher than is typical for two-photon 
interactions. Particularly at the LHC, where M/7 is very small, pt should be an effective 
tool for separating the two event classes. 



VI. CONCLUSIONS 



We have calculated the total cross sections and rapidity and transverse momentum distri- 
butions for two-photon production of various final states with and without nuclear breakup. 
The cross-sections for 77 interactions accompanied by XnXn and lnln dissociation are 
about 1/10 and 1/100 of that for the unaccompanied 77 interactions. The nuclear breakup 
tagged events have different rapidity and px distributions from the un-tagged events, and 
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can be used to explore the effects of different photon spectra and b distributions. The typi- 
cal pt of 77 final states is several times smaller than for comparable coherent photonuclear 
interactions, so pr cuts may be an effective tool for separating the two interactions. 

This work was supported by the U.S. Department of Energy under Contracts No. DE- 
AC-03076SF00098 and DE-AC02-98CH10886. 
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